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Abstract 



k> i In order to obtain polarised parton densities we have made next to 

^ '_ leading order QCD fit using experimental data on deep inelastic struc- 

ture functions on nucleons. This fit is compared with the updated fit 
to corresponding spin asymmetries. We get very similar results for 
all fits also for different data samples. It seems that only polarised 
parton densities for non-strange quarks An and Ad are relatively well 
determined from the present polarised deep inelastic experiments. In- 
tegrated gluon contribution at Q'^ = 1 GeV^ is, as in our previous fits, 
very small. 



In order to get polarised parton (i.e. quark and gluon) densities one uses 
data on deep inelastic scattering of polarised lepton on polarised nucleon 
targets. Quite a lot of data exist for such scattering on different nucleon 
targets. The data come from experiments made at SLAC [1-10], CERN 



[11-16] and DESY [0, |T8l. Recently the data from SLAC from E155 |0 
experiment on protons has been published. The experimental groups present 
data for spin asymmetries as well as for polarised structure functions. 



The analysis of the EMC group results [|TT| started an interest in studying 



such data. So called "spin crisis" was connected with the fact that only very 
little of spin of nucleon was carried by quarks. The suggestion from ref. |]19 



was that polarised gluons may be responsible for this phenomena. Many next 
to leading (NLO) order QCD analysis [20-27] were performed and polarised 
parton distributions were determined. The main purpose of this paper is 
to use data for polarised structure functions 5'i(x, Q^) on proton, neutron 
and deuteron targets in order to determine polarised parton distributions. 
This fit will be compared with our updated (in which we take into account 
recently published data on protons from E155 |T0[ experiment in SLAC) 



fits to spin asymmetries. It was advocated by us [^ and [^ that using 
spin asymmetries for determination of polarised parton densities one avoids 
the problem with higher twist contributions. On the other hand it is the 
polarised structure fuctions and polarised parton distributions that we want 
to determine. 

We compare these two ways of making fits with similar technical as- 
sumptions and the same parton functions used for fitting. We will see that 
both methods give very similar results for parton distributions. Our pre- 
vious p6| , p7| important conclusion that the integrated gluon contribution 
is negligible a.t Q^ = 1 GeV^ does not change. Most of the groups used 
experimental data for spin asymmetries to determine polarised parton dis- 
tributions. In addition to spin asymmetries one has also experimental data 
for polarised structure functions calculated from the spin asymmetries in a 
specific way chosen by an experimental group. There were also several fits 
to the data on polarised structure functions [^ ^. As in |]2B|, ^ we will 



make fits to two samples of the data. In the first group we will have data 
for the same x (strictly speaking for the near values) and different Q^ and in 
the second the "averaged" data where one averages over Q^ (the errors are 
smaller and Q^ dependence is smeared out). In most of the fits to experi- 
mental data only second sample (namely with averaged Q^ dependence) was 



used. Our fits use the both sets of data. The data for polarised structure 
functions are usually given only for averaged sample of data (the exception 
is E155 experiment for deuteron and proton data). 

Experiments on unpolarised targets provide information on the unpo- 
larised quark densities q{x, Q"^) and G(a;, Q"^) inside the nucleon. These den- 
sities can be expressed in term of g^(x,Q^) and G^{x,Q'^), i.e. densities of 
quarks and gluons with hehcity along or opposite to the helicity of the parent 
nucleon. 

g = g+ + g-, G = G+ + G' . (1) 

q stands for quark and antiquark contributions. 

The polarised parton densities, i.e. the differences of g"*", q~ and G^ , G~ 
are given by: 

Ag = g+-g~, /\G = G+ - G' . (2) 

We will try to determine q^{x, Q^) and G^{x, Q"^), in other words, we will 
try to connect unpolarised and polarised data. 

In our fits we will use functions for the polarised parton densities that are 
suggested by the fit to unpolarised data. We risk that asymptotic behaviour 
of our parton distributions is not correct but it seems to us not so important 
when we limit ourselves to the measured region of x. Maybe it is also not 
bad to use parametrisation completely different from other groups to check 
how the results depend on that. Our starting point are the formulas for 
unpolarised quark and gluon distributions gotten (at Q^ = 1 GeV^) from 
the fit performed by Martin, Roberts, Stirling and Thorne |^ (they use 
A^^ = 0.3 GeV and a.(M|) = 0.120) 

We will not use small and high x behaviour of unpolarised parton dis- 
tributions as fitted parameters as some other groups do. We will split q 
and G, as was already discussed in ref. P^, P^, into two parts in such a 



manner that the distributions q^{x, Q"^) and G'^{x, Q^) remain positive. Our 
polarised densities for quarks and gluons are parametrised as follows: 



^u^{x) = x-^-^^^\l - xf^^^ai + a2v^ + a^x), 
Ad,{x) = x~°-'"^(l - xf-^'^bi + b2V^ + hx), 

2Au{x) = OAAM{x) - A6{x), (3) 

2Ad{x) = 0.4AM(x) + AS{x), 



2As(x) = 0.2AM,(x), 

AG{X) = X-°-°^29^^ _ a;)6.587(^^ ^ t/2V^+ t/gx). 

whereas for the antiquarks ( and sea quarks ): 

AM(x) = x-'-^'^^il - xy-''^ici + C2v^), 
AM, = x-'-'''\l - x)^-«o«(ci, + C2.v^), (4) 

A(5(x) = x°-^^=^(l - x)9-«°^C3(l + 9.987X - 33.34x2), 

We also have 



Am = Auy + 2Au, 

Ad = Ad^ + 2 Ad, (5) 

As = 2As, 



We use as before [^, ^ additional independent parameters for the strange 
sea contribution with the same as for non-strange sea functional dependence. 
Maybe not all parameters are important in the fit and it could happen that 
some of the coefficients in eq.(3,4) taken as free parameters in the fit are 
small or in some sense superfiuous. Putting them to zero (or eliminating 
them) increase x^ only a little but makes the quantity x^/^df smaller. We 
will see that that is the case with some parameters introduced in eq. (3,4). 
In order to get the unknown parameters in the expressions for polarised 
quark and gluon distributions (eqs.(3,4)) we calculate the spin asymmetries 
(starting from initial Q^ = 1 GeV^) for measured values of Q^ and make a 
fit to the experimental data on spin asymmetries for proton, neutron and 
deuteron targets. The spin asymmetry Ai{x,Q^) can be expressed via the 
polarised structure function gi{x,Q'^) as 

where R = [-^2(1 +7^) ~ 2x-Fi]/2xFi whereas Fi and F2 are the unpolarised 
structure functions and 7 = 2Mx/Q (M stands for proton mass). We will 
take the new determined value of R from the [^]. The factor (1 + 7^) plays 
non negligible role for x and Q^ values measured in SLAG experiments. In 



calculating gi{x, Q"^) and F2{x, Q"^) in the next to leading order we use proce- 



dure described in [^] following the method described in ||20|, ^ performing 
calculations with Mellin transforms and then calculating Mellin inverse. Hav- 
ing calculated the asymmetries according to equation (6) for the value of Q"^ 
obtained in experiments we can make a fit to asymmetries on proton, neu- 
tron and deuteron targets. The other possibility is to use directly the data 
for the polarised structure function f7i(x, Q^) (the values of gi were given for 
the averaged values of Q"^) on proton, neutron and deuteron targets to deter- 
mine unknown coefficients in expressions for polarised parton distributions. 
The problem is that different experimental groups used their own specific 
methods to obtain the values of polarised structure functions gi{x, Q"^) from 
the measured asymmetries. At the end we want to know polarised struc- 
ture functions and polarised parton distributions. In order to calculate them 
from spin asymmetries we have to choose what shall we take for Fi{x, Q"^) or 
F2{x,Q^) and R{x,Q^). As was already mentioned we calculated F2{x,Q'^) 
in NLO for actual values of x and Q^ using quark and gluon contributions 
for Q^ = 1 GeV^ given by MRST [^ . The values of R were in the earlier 
fits taken from Whitlow ||32| and later from E143 group |0 .We have treated 



all experiments in the same way. There is also a problem of higher twist 
corrections (power low corrections to R were included). We will not include 
the higher twist corrections because of still big experimental errors. The 
spread of the results could be a measure of uncertainties in both methods. 
We will compare fits using determination of parameters from polarised struc- 
ture functions and from spin asymmetries. As will be seen later the results 
obtained by two methods are very similar. 

It was already seen from our previous paper p^ that gi{x, Q"^) calculated 



from spin asymmetries fits not bad the data points for gi. Data points for 
polarised structure functions were given for averaged data set so it is natural 
to compare fit to gi (called by us fit g) with the fit to spin asymmetries (fit Ai) 
with the same number of points. We will take into account in this case 197 
data points (We will take E155 proton and deuteron data without averaging). 
These fits will be compared with the fit (fit A2) to spin asymmetries for non 
averaged data where we take into account 465 experimental data points. At 
the beginning we will not put any constrains from hyperon decays. Later we 
will also not fix og = Am + Ad — 2As value but we will add experimental 
point Og = 0.58 ±0.1 with enhanced (to 3a) error. That means we will 
simply add to x^ corresponding to experimental points for spin asymmetries 



the term connected with experimental point from hyperon decays. We will 
discuss how this additional experimental point infiuences our results. 

Not all parameters are important in the fits. It seems that some of the 
parameters of the most singular terms are superfiuous and we can eliminate 
them. We will put di = (such assumption gives that 5G/G ~ x^^"^ for 
small x), bi = (the most singular term in Ad^) and assume Cis = Ci (i.e. 
the most singular terms for strange and non-strange sea contributions are 
equal). Fixing these four parameters in the fit practically does not change 
the value of x^ but improves x^/^df- We also have to make some remarks 
about parameters C3 and ^2- In the fit for gi these parameters are important 
i.e. when we eliminate them x^ per degree of freedom increases. That is 
not the case in fits for spin asymmetries. In this case these parameters are 
superficial (that means for example that splitting u — d is not well defined 
by data to spin asymmetries. We will leave this parameters for comparison 
with fit to (?i (do not eliminate them) but they could be not well determined 
and cause some artificial shifts in other parameters. 



Table 1. The parameters of our three fits calculated at Q^ 
together with x^ P^i" degree of freedom 



lGeV2 



fit 


ai 


^2 


04 


&2 


h 


Cl 


&tg 


0.61 
±0.0 


-7.05 
±0.23 


17.1 
±0.23 


-2.02 
±0.0 


0.34 

±0.24 


-0.34 
±0.03 


fit Ai 


0.56 
±0.13 


-5.50 

±1.22 


14.7 
±1.66 


-1.67 
±0.03 


-0.173 
±0.26 


-0.338 
±0.10 


fit A2 


0.49 
±0.13 


-5.34 
±1.20 


14.66 
±1.63 


-2.02 
±0.0 


0.35 

±0.25 


-0.32 
±0.10 



fit 


C2 


C2s 


C3 


d2 


^3 


x'INdf 


fit(? 


4.15 
±0.0 


4.15 
±0.0 


-1.05 
±0.56 


-29.0 

±8.6 


87.1 
±36.1 


0.87 


fit Ai 


3.23 

±0.80 


4.15 
±0.22 


-0.617 
±0.58 


-15.4 
±0.22 


42.2 
±15.0 


0.81 


fit A2 


3.69 

±0.79 


4.15 
±0.20 


-0.39 

±0.48 


-14.0 
±0.04 


27.0 
±11.4 


0.84 



In the Table 1 we present the values of parameter from the fit to the data 
on polarised structure functions and spin asymmetries for averaged and non 
averaged data together with x^/^df values. 

For example in the case of first row in table 1 corresponding to the fit 
to polarised structure functions the obtained quark and gluon distributions 
lead for {Q^ =1 GeV^) to the following integrated (over x) quantities: Au = 
0.72, Arf = -0.64, As = 0.05, Am^ = 0.54, Arf^ = -0.65, 2Am = 0.18, 2AJ = 
0.01. 

We have positively polarised sea for up and down quarks and positively 
polarised sea for strange quarks. The gluon polarisation is small. The value of 
03 = 1.36 was not assumed as an input in the fit (as is the case in nearly all fits 
[ P^ ) and comes out slightly higher than the experimental value. The value 
of Og = —0.01 is completely different from the experimental figure. Taking 
into account that fits to polarised structure functions and spin asymmetries 
use different methods to calculate F2{x, Q"^) and R{x, Q"^) there is no reason 
to expect that they give exactly the same results. The obtained values of 
parameters are very close and practically within experimental errors. The 
parameters calculated in fits to spin asymmetries (line 2 and 3) are closer 
in comparison with fit to gi (line 1) but are not identical. The spread of 
parameters measures small differences in the Q^ evolution, differences in 
experimental errors and influence of our specific functions used in fits. 



As was already mentioned in p6| the asymptotic behaviour at small x of 
our polarised quark distributions is determined by the unpolarised ones and 
hence do not have the expected theoretically Regge type behaviour. Some 
of the quantities specially integrated sea contributions and also some valence 
contributions in our fit change rapidly for x < 0.003. that is not something 
that we expect from Regge behaviour with small exponent. 

Hence, we will present quantities integrated over the region from a;=0.003 
to x=l (it is practically integration over the region which is covered by the 
experimental data, except of non controversial extrapolation for highest x). 
The values of integrated quantities in the measured region we consider as 
more reliable then those in the whole region. 

The corresponding quantities for three our fits together with x^/Ndf are 
presented in table la. 
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Figure 1: The quark and gluon densities for up and down flavour versus x 
gotten from the fit g (solid line) and Ai (dashed line). 



Table la. The values of quark and gluon polarisations at Q^ 
for our three fits 



IGeV^ 



fit 


An 


Ad 


Au^ 


Ad^ 


2Au 


2Ad 


2As 


AG 


fit (7 


0.74 


-0.49 


0.44 


-0.63 


0.30 


0.14 


0.11 


0.15 


fit Ai 


0.75 


-0.48 


0.57 


-0.57 


0.18 


0.09 


0.11 


0.01 


fit A2 


0.76 


-0.47 


0.54 


-0.63 


0.22 


0.16 


0.12 


-0.19 



From the table we see that there are changes in valence and sea contribu- 
tions in different fits but the values of Am and Ad practically do not differ. 
We use the parametrisation where the most singular term in sea contribu- 
tion is very similar to valence quark terms and that maybe is the reason why 
splitting into valence and sea contribution is fragile and changes for different 
fits but in the case of Au and Ad do not differ much. In the first fit at 
Q2 = I GeV^ we get F? = 0.124, F^ = -0.052, as = 1.23 and AS = 0.36 
comparing with the second fit to averaged spin asymmetries where we have 
F? = 0.125, F^ = -0.051, ag = 1.24 and AS = 0.38. These results are very 
close. The value of 03 in the measured region without any assumption comes 
out close to the value measured in hyperon decays. 

For illustration we present in Fig.l the distributions Am and Ad for our 
sets of parameters calculated from polarised structure functions and the sam- 
ple of averaged spin asymmetries data. The corresponding values for Au^ 
and Ady differ much stronger. In our previous paper ||2^ we already pre- 
sented how the values of polarised structure functions for proton, deuteron 
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Figure 2: The comparison of our predictions for gi{x,Q'^) versus x with 
experimental data from different experiments. Solid curve is gotten from fit 
g, the dashed one is calculated using the parameters of fit A\. 



and neutron calculated from the fits to spin asymmetries compare with the 
experimental data. To see what is the difference in the values fitted directly 
to the polarised structure functions and the values calculated from the fit 
to spin asymmetries we present in Fig. 2 the corresponding curves in com- 
parison with experimental points for gf, gf and (/" at the values of Q^ in 
corresponding experiment. As an example we show comparison with experi- 
mental points for polarised structure functions (^f, gf from SMC from CERN 
and E143 from SLAC and (7" from E142 and E154 from SLAC. There are 
some differences but they are not big in comparison with experimental errors. 
The comparisons for other experimental sets look very similar. It is of course 
not astonishing that the fitted curves with the parameters from the first fit 
are closer to experimental values. 

As we already pointed out before we have not made any assumptions 
about Og. We got from the fits that the value of ag is near zero very far from 
the experimental value and we got positive values for As. The value of 03 
that also was not constrained in the fit is close to experimental value (in the 
measured region of x). In order to make more direct comparison with other 
fits as before we will also not fix ag = A-u + Arf — 2As value but we will 
add experimental point ag = 0.58 ±0.1 with enhanced (to 3o") error. That 
means we will simply add to x^ corresponding to experimental points for 
spin asymmetries the term connected with experimental point from hyperon 
decays. The parameters of our three new fits (called g' , A[ and A2 ) are now 
presented in Table 2 and results in the measured region in 0.003 < x < 1 in 
the Table 2a. 
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Table 2. The parameters of three new fits calculated at Q^ = 1 GeV^ 



fit 


ai 


02 


a4 


&2 


&3 


Cl 


Rig' 


0.61 


-6.84 


16.83 


-1.86 


0.13 


-0.31 


fit A[ 


0.56 


-5.51 


14.73 


-1.65 


-0.21 


-0.34 


fit A'2 


0.50 


-5.39 


14.72 


-1.98 


0.29 


-0.33 



fit 


C2 


C2s 


C3 


^2 


d. 


xVNdf 


atg' 


4.15 


-0.54 


-1.04 


-34.5 


102.2 


0.88 


fit A[ 


3.67 


-0.63 


-0.64 


-15.8 


43.5 


0.80 


fit A'^ 


4.15 


-0.56 


-0.44 


-14.2 


27.6 


0.84 



There are some small changes in the parameters in comparison to the 
fits g, Ai, A2 and the biggest change is in C2s the parameter responsible for 
the strange sea. Strange quark contribution is not well determined by the 
polarised deep inelastic data alone and it is easy by additional experimental 
point on as from hyperon decays to shift the value of as from nearly zero 
to correct experimental value with only small changes in non strange parton 
parameters. Comparing Table la and Table 2a we can see what is the infiu- 
ence of this additional experimental point as on integrated parton densities 
for our three fits. This additional experimental point causes shifts of inte- 
grated parton values. As is seen from Table 2a the valence non-strange quark 
distributions nearly cancel and the value of ag = A-u + Ad — 2As = 0.58 is 
built up from relatively high sea contributions. The values of Ti, F" and as 
do not change in comparison to previous fits. One can easily calculate from 
the Table 2a that AS = 0.24 in the fits g' and A[. 



Table 2a. The values of quark and gluon polarisations at Q'^ = 1 GeV^ for 
three new fits (where one includes experimental point from hyperon decays). 



fit 


Am 


Ad 


Au,, 


Ad, 


2Au 


2Ad 


2As 


AG 


mg' 


0.79 


-0.44 


0.47 


-0.60 


0.32 


0.16 


-0.11 


0.16 


fit a; 


0.80 


-0.44 


0.57 


-0.57 


0.23 


0.13 


-0.12 


0.01 


fit A'2 


0.80 


-0.43 


0.54 


-0.62 


0.26 


0.19 


-0.11 


-0.19 
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However it is not clear whether the general conclusion that the sea con- 
tributions for quarks ( both non-strange and strange) are very big in the 
measured x region is correct. It is specific for our model that the leading 
singularity for polarised quark valence and sea contributions are comparable. 
That means that splitting into valence and sea contributions could be not 
well determined in our fits. That of course could be connected with the func- 
tional form of polarised parton densities used by us but what is important 
is to reproduce the results of the experiment in the measured region. That 
means that only Am and Ad values can be determined using our parametri- 
sation from the polarised deep inelastic data and not valence and sea con- 
tributions separately. The value of As is determined only with additional 
Og value from hyperon decays. In our model AG is small. The integrated 
values Am = 0.80, Ad = —0.44 are actually not unexpected. Similar values 
follow from the other models. In other fits I^D], ^ with completely different 
assumptions for example when we assume the values of 03 and Og by fixing 
the parameters of the fitted parton distributions (normalisation constants) 
and with the assumption of SU{3) symmetry for quark sea we get (using 
completely different parametrisation from that we use for polarised parton 
densities) : 

Auy - A4 = 1.26, 

Auy + Ad^ = 0.58 (7) 

If follows that Auy = 0.92, A4 = -0.34 and in order to get AE = 0.20 
[m we get 2Au = 2Ad = 2As = -0.13 and that means Am = 0.79, 6d = 
—0.47 ( the values not very different from our values). In such models sea 
contribution is relatively big and negative contrary to our model where we 
have at least in the measured region big and positive sea contribution. This 
type of splitting into big valence and relatively big negative sea contribution is 
caused by the assumptions of the model. We have not made such assumptions 
taking ag as additional experimental point. Our solution with relatively small 
valence contribution and relatively large positive sea contributions is different 
but the values of Am and Ad in both models are very close. It seems that our 
assumptions are less restrictive. The fact that we can get completely different 
splitting into valence and sea contributions using the same experimental data 
shows that this splitting is not well determined by experimental data what 
is more reliable are distributions of Am and Ad and their integrated values. 

12 



We have made fits to data on polarised structure functions on proton, 
neutron and deuteron targets and we liave determined polarised parton dis- 
tributions. These fits are compared with our previous fits for corresponding 
asymmetries (improved by usage of recent data from El 55 proton experi- 
ment in SLAC). As a check we also have made fits to spin asymmetries for 
all (non averaged in Q'^) data on spin asymmetries. These fits lead to very 
similar results with small integrated gluon contribution. The fits were made 
without inclusion of information on 03 and ag from hyperon decays and then 
repeated with additional experimental point on ag- In the first case ag is close 
to zero and As is positive. In the second case additional experimental point 
on as changes practically only parameters of strange quark and causes small 
shifts in other parameters. The value of a^ at least in the measured region 
of X without any assumptions comes out very close to experimental value. 
It seems that with the parametrisation used by us only Au{x) and Ad{x) 
are well determined (not the splitting into valence and sea parts). Polarised 
strange quark distributions, gluon distributions and also u — d splitting are 
not well determined by polarised deep inelastic experimental data. 
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